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ABSTRACT 

We report extensive spectroscopic and differential V-band photometric observations of the 18.4-day 
detached double-lined eclipsing binary LV Her (F9 V), which has the highest eccentricity (e ~ 0.613) 
among the systems with well-measured properties. We determine the absolute masses and radii of 
the components to be M 1 = 1-193 ± 0.010 M Q , M 2 = 1.1698 ± 0.0081 M Q , R x = 1.358 ± 0.012 i? Q , 
and i?2 = 1-313 ± 0.011 Rq, with fractional errors of 0.9% or better. The effective temperatures are 
6060 ± 150 K and 6030 ± 150 K, respectively, and the overall metallicity is estimated to be [m/H] 
= +0.08 ± 0.21. A comparison with current stellar evolution models for this composition indicates 
an excellent fit for an age between 3.8 and 4.2 Gyr, with both stars being near the middle of their 
main-sequence lifetimes. Full integration of the equations for tidal evolution is consistent with the 
high eccentricity, and suggests the stars' spin axes are aligned with the orbital axis, and that their 
rotations should be pseudo-synchronized. The latter prediction is not quite in agreement with the 
measured projected rotational velocities. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: evolution — stars: funda- 
mental parameters — stars: individual (LV Her) 



1. INTRODUCTION 

The photometric variability of the eclipsing binary 
star LV Her (TYC 2076-1042-1; a = 17 h 35 m 32*40, 
5 = +23° 1 0' 30"6, J200 0; Sp T F9, V = 11.02) was an- 
no unced bvlHo ffmeister (|1935l ) . Its period w as estimated 
bv | Zessewitschl (|1944l ) as 2.634 days , by iZessewitschl 
(119541) as 5.2" 674 days, bv lPopperl (fl99l as 9.218 days, by 
iTorresI (pOOOl ) as 18.1312 days, and by lTorres et alJ (pOOlh 
as 18.4359350 days. All photometric and spectroscopic 
work performed since this last study has confirmed that 
it is basically correct. The prior erroneous period values 
resulted from the very long orbital period, the similar 
eclipse depths (0.68 mag and 0.66 mag in V; see below), 
the extremely high eccentricity of the orbit (0.613), and 
the resultant shift of secondary eclipse to a phase of 0.86. 
This sort of erroneous p eriod estimation is not un usual 
in this kind of situation (|Lacv et al.ll2004allb1 . [200l) . 

Largely because of this uncertainty, no determination 
of the absolute dimensions of the system has been made 
until now. In addition to an accurate determination of 
the eclipse ephemeris in § 12. 2[ this work presents new 
high-quality photometric and spectroscopic observations 
(§ 12.11 and § 12. 3| that yield masses and radii on a par 
with the best determinations for eclipsing binaries to 
date (§HJ). In § !5.1l we compare these determinations with 
models of stellar structure and stellar evolution. The 
high eccentricity and long period of LV Her make it an 
interesting case for comparison with tidal theory, which 
we present in § 15.21 We discuss our results in §[6j 



2. OBSERVATIONS AND REDUCTIONS 
2.1. Differential and absolute photometry 
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Our photometric work on LV Her began in 2001 and 
is based on observations obtained with two robotic in- 
struments: the U RSA telescope at Kimpel Observatory 
(jLacv et all l2001h on the campus of the University of 
Arkansas, and a robotic telesc ope at the NF Observa- 
tory (NFO: iGrauer et al.lf2008h near Silver City NM. 

Kimpel Observatory (ursa.uark.edu) consists of a 
Meade 10-inch f/6.3 LX-200 telescope with a Santa Bar- 
bara Instruments Group ST8 CCD camera (binned 2x2 
to produce 765 x 510 pixel images with 2.3 arcsec square 
pixels) inside a Technical Innovations Robo-Dome, and 
controlled automatically by an Apple Macintosh G4 com- 
puter. The observatory is located on top of Kimpel Hall 
on the Fayetteville campus of the University of Arkansas, 
with the control room directly beneath the observatory 
inside the building. Exposures of 60 or 120 seconds made 
through a Bessell V filter (2.0 mm of GG 495 and 3.0 mm 
of BG 39) were read out and downloaded by ImageGrab- 
ber (camera control software written by J. Sabby) to the 
control computer over a 30 second interval, then the next 
exposure was begun. The observing cadence was there- 
fore about 90 to 150 sec per observation. The variable 
star would sometimes be monitored continuously for 4-6 
hours. LV Her was observed by URSA on 156 nights dur- 
ing parts of nine observing seasons from 2001 February 18 
to 2009 April 8. The total number of URSA observations 
is 6690. The frames were analyzed by a virtual measur- 
ing engine application written by Lacy that flat-fielded 
the images, automatically located the variable and com- 
parison stars in the image, measured their brightnesses, 
subtracted the corresponding sky brightness, and cor- 
rected for the differences in airmass between the stars. 
Extinction coefficients were determined nightly from the 
comparison star measurements. They averaged 0.25 
mag/airmass. LV Her is TYC 2076-1042-1. The com- 
parison stars were TYC 2076-0580 -1 (comp, V T = 11.12, 
as listed in the Tycho-2 Catalogue: iHpg et al. 2000^, and 
TYC 2076-1387-1 (ck, V T = 11.53). Both comparison 
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TABLE 1 

URSA DIFFERENTIAL V-BAND MAGNITUDES OF 

LV Her. 



HJD-2,440,000 


AV 


Orbital phase 


51958.87391 . . . 


-0.195 


0.15136 


51958.87498 . . . 


-0.147 


0.15141 


51958.87603 . . . 


-0.174 


0.15147 


51958.87708 . . . 


-0.179 


0.15153 


51958.87812... 


-0.183 


0.15158 



Note. — Table [T] is available in its entirety 
in the electronic edition of the Astronomical 
Journal. A portion is shown here for guidance 
regarding its form and contents. 



TABLE 2 

NFO DIFFERENTIAL V-BAND MAGNITUDES 

of LV Her. 



HJD-2,440,000 


AV 


Orbital phase 


53430.93043 . . . 


1.003 


0.99840 


53430.93323 . . . 


1.010 


0.99855 


53430.93602 . . . 


1.035 


0.99870 


53430.93878 . . . 


1.055 


0.99885 


53430.94158 . . . 


1.068 


0.99900 



Note. — Table [2] is available in its en- 
tirety in the electronic edition of the Astro- 
nomical Journal. A portion is shown here 
for guidance regarding its form and con- 
tents. 

stars are within 10 arcmin of the variable star (var). 
The comparison star magnitude differences (comp— ck) 
were constant at the level of 0.013 mag (standard devi- 
ation within a night), and 0.010 mag for the standard 
deviation of the nightly mean magnitude difference. The 
URSA differential magnitude (var— comp) of the variable 
star was referenced only to the magnitude of the com- 
parison star, comp. The resulting 6690 V magnitude dif- 
ferences (var— comp) are listed in Table [T] (without any 
nightly corrections) and plotted in Figures [T] [3] (after the 
nightly corrections discussed below have been applied). 
The eclipse depths are 0.68 mag for the primary and 
0.66 mag for the secondary. The precision of the vari- 
able star differential magnitudes is about 0.016 mag. 

The other telescope we used is the NFO WebScope, a 
refurbished 24-inch Group 128 Cassegrain reflector with 
a 2K x 2K Kodak CCD camera, located near Silver City, 
NM. Observations consisted of 120 second exposures also 
through a Bessell V filter. LV Her was observed by NFO 
on 123 nights during parts of five observing seasons from 
2005 March 1 to 2009 April 24, yielding 2946 observa- 
tions. Extinction coefficients were determined nightly 
from the comparison star measurements. They averaged 
0.18 mag/airmass. The same comparison stars were used 
as those of the URSA telescope. The comparison star 
magnitude differences (comp— ck) were constant at the 
level of 0.007 mag (standard deviation within a night), 
and 0.020 mag for the standard deviation of the nightly 
mean magnitude difference. The differential magnitude 
(var— comps) of the variable star was referenced to the 
magnitude corresponding to the sum of the intensities of 
the comparison star, comp, and the check star, ck. The 
resulting 2946 V magnitude differences (var— comps) are 
listed in Table [2] (without nightly corrections) and plot- 



ted in Figure HH3] (including nightly corrections; see be- 
low) . The precision of the variable star differential mag- 
nitudes is about 0.010 mag. We noticed early on dur- 
ing the observations that the NFO magnitudes showed a 
small but significant offset from night to night, on the or- 
der of a hundredth of a magnitude. The origin of the off- 
set is a variation in responsivity across the field of view of 
the NFO combined with imprecise centering from night 
to night. These variations are a well-known effect of the 
optics when using wide-field imaging telescopes such as 
the NFO. We have removed most of this variation by 
using dithered exposures of open star clusters to mea- 
sure this variation, fitting a 2-D polynomial (see lSelmanl 
2004), and subtracting the variation during initial re- 
ductions (photometric flat). The URSA observations, on 
the other hand, show this kind of effect to a very much 
smaller extent. We have removed these nightly offsets 
before further analysis by using a procedure discussed 
below in §[3] 

Absolute photometry of LV Her is available in the lit- 
erature from several sources, and in several photomet- 
ric systems including Stromgren, 2MASS, Tycho-2, and 
Johnson-Cousins. We collect these measurements in Ta- 
ble [H Color indices formed from these magnitudes can 
be used to estimate a mean effective temperature for the 
binary, which we discuss below in §2J The interstellar 
reddening can be estimated, for examp le, by comparison 
with t he standard Stromgren indices of lPerrv fe; Johnson! 
(1982), giving E{b-y) = 0.007 mag, or E(B-V) « 0.01. 
Another estimate is available from the reddening maps o f 
the sky near the variable star from lSchlegel et all (1998), 
which results in a larger value of E{B — V) = 0.057 mag 
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Fig. 1. — URSA and NFO V-band photometry for LV Her, shown 
with our best fit model described in §[3] Residuals are shown at 
the bottom, with NFO displaced for clarity. 
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Fig. 2.— Enlarged view of the URSA and NFO V-band light 
curves for LV Her near the primary eclipse, shown with our best fit 
model described in §|3] Residuals are shown at the bottom, with 
NFO displaced for clarity. 



TABLE 3 

Absolute photometry for LV Her (combined 
light). 



Passband Value Reference 



y 11.02 ± 0.02 a 1 

u- b 1.552 ± 0.030 a 1 

v-b 0.571 ± 0.020 a 1 

b-y 0.367 ± 0.015 a 1 

V 11.055 ± 0.040 2 

I c 10.367 ± 0.058 2 

J 9.905 ± 0.018 a 3 

H 9.665 ± 0.025 a 3 

K s 9.631 ± 0.018 a 3 

B T 11.832 ± 0.066 4 

V T 11.045 ± 0.055 4 



Note. — References: 1 -IHilditch k. H ill (1975); 



2 - TASS IDroege et al., 20071) ; 3 - 2MASS; 4 - 
Tycho-2 IIH0g et alj|200ch . 

a Single measurement obtained out of eclipse. 

at the estimated distance of the eclipsing binary. We 
adopt for further use the straight average of these two 
reddening estimates, E(B — V) — 0.03 ± 0.03 mag, with 
a conservative error. 
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Fig. 3.— Enlarged view of the URSA and NFO V-band light 
curves for LV Her near the secondary eclipse, shown with our best 
fit model described in §[3] Residuals are shown at the bottom, with 
NFO displaced for clarity. 

light (37 for the primary, 30 for the secondary) are listed 
in Table 21 All available minima have been used to esti- 
mate an eclipse ephemeris by a least squares technique. 
When available, the uncertainties have been adopted as 
published. For the others (visual and photographic mea- 
surements), uncertainties have been estimated by iter- 
ations requiring the reduced x 2 to be near unity, sep- 
arately for each class of measurement technique. For 
the visual meaurements we obtain a — 0.017 days, and 
for the photographic timings a = 0.053 days. Separate 
solutions with the primary and secondary data give pe- 
riods that are not significantly different. For the final 
ephemeris we enforced a common period, and obtained: 

Min 1 = 2,453,652.19147(10)+ 18.4359535(19)£ 
Min 11 = 2,454,165.86045(12)+ 18.4359535(19)£ 

where the reference epochs were chosen so as to min- 
imize the correlations with the period and with each 
other. This ephemeris is the basis of the orbital phases we 
cite below. The phase of secondary eclipse is 0.86235 ± 
0.00001. A plot of the O-C diagram from this fit is 
shown in Figure [H and the residuals are included in Ta- 
ble [4] It is quite remarkable how much the accuracy has 
improved over the last century. 



2.2. Ephemeris 

Since the work of lTorres et al.l (|2001[ ). additional times 
of minimum for LV Her have been reported in the litera- 
ture, and some recent ones from our WebScopes are avail- 
able that have not yet been published. The previously 
published ones and these additional times of minimum 



2.3. Spectroscopy 

LV Her was placed on the observing list at the Harvard- 
Smithsonian Center for Astrophysics (CfA) on 1991 May 
31st, and was monitored spectroscopically for exactly 
11 years with an echelle instrument on the 1.5m Till- 
inghast reflector at the F. L. Whipple Observatory on 



Torres et al. 



TABLE 4 
Times of eclipse for LV Her. 



HJD-2,400,000 


Year 


Type 


Instrument 


a (d) 


O-C (d) 


Reference 


15929.712 


1902.4907 


2 


Pg 


0.053 


+0.01915 


1 


16224.798 


1903.2986 


2 


Pg 


0.053 


+0.12989 


1 


16669.658 


1904.5165 


1 


Pg 


0.053 


-0.01072 


1 


17351.799 


1906.3841 


1 


Pg 


0.053 


-0.00000 


1 


18605.473 


1909.8165 


1 


Pg 


0.053 


+0.02916 


1 


23601.681 


1923.4954 


1 


Pg 


0.053 


+0.09376 


1 


24062.505 


1924.7570 


1 


Pg 


0.053 


+0.01892 


1 


26032.643 


1930.1510 


2 


Pg 


0.053 


+0.04762 


1 


26032.644 


1930.1510 


2 


Pg 


0.053 


+0.04862 


1 


26090.482 


1930.3093 


1 


Pg 


0.053 


+0.04103 


1 


26219.577 


1930.6628 


1 


Pg 


0.053 


+0.08436 


1 


26901.699 


1932.5303 


1 


Pg 


0.053 


+0.07608 


1 


27212.514 


1933.3813 


2 


Pg 


0.053 


+0.01760 


1 


27636.462 


1934.5420 


2 


Pg 


0.053 


-0.06134 


1 


27657.429 


1934.5994 


1 


Pg 


0.053 


-0.06802 


1 


28281.830 


1936.3089 


2 


Pg 


0.053 


+0.04829 


1 


28429.2740 .... 


1936.7126 


2 


vis 


0.017 


+0.00466 


2 


28487.1220 .... 


1936.8710 


1 


vis 


0.017 


+0.00708 


2 


30254.522 


1941.7098 


2 


Pg 


0.053 


+0.09326 


1 


31268.42 


1944.4857 


2 


vis 


0.017 


+0.01382 


3 


31289.35 


1944.5431 


1 


vis 


0.04 


-0.02986 


3 


31326.29 


1944.6442 


1 


vis 


0.017 


+0.03823 


3 


31342.16 


1944.6876 


2 


vis 


0.017 


+0.01001 


3 


34626.332 


1953.6792 


1 


Pg 


0.053 


+0.04456 


1 


37444.539 


1961.3950 


2 


Pg 


0.053 


+0.08839 


1 


38147.573 


1963.3198 


1 


Pg 


0.053 


+0.01843 


1 


38587.449 


1964.5242 


2 


Pg 


0.053 


-0.03072 


1 


40744.498 


1970.4298 


2 


Pg 


0.053 


+0.01171 


1 


41060.502 


1971.2950 


1 


Pg 


0.053 


+0.06678 


1 


41982.251 


1973.8186 


1 


Pg 


0.053 


+0.01810 


1 


42664.347 


1975.6861 


1 


Pg 


0.053 


-0.01618 


1 


46001.268 


1984.8221 


1 


Pg 


0.053 


-0.00276 


1 


46941.483 


1987.3963 


1 


Pg 


0.053 


-0.02139 


1 


47381.384 


1988.6006 


2 


Pg 


0.053 


-0.04555 


1 


48100.425 


1990.5693 


2 


Pg 


0.053 


-0.00674 


1 


48487.60 


1991.6293 


2 


vis 


0.017 


+0.01328 


1 


48508.60 


1991.6868 


1 


vis 


0.04 


+0.03956 


1 


48545.52 


1991.7879 


1 


vis 


0.04 


+0.08765 


1 


48745.70 


1992.3359 


2 


vis 


0.017 


+0.00989 


1 


48948.50 


1992.8912 


2 


vis 


0.017 


+0.01440 


1 


49098.540 


1993.3020 


1 


Pg 


0.053 


+0.02904 


1 


49206.61 


1993.5978 


2 


vis 


0.017 


+0.02105 


1 


49667.49 


1994.8597 


2 


vis 


0.017 


+0.00221 


1 


49688.48 


1994.9171 


1 


vis 


0.017 


+0.01853 


1 


49925.65 


1995.5665 


2 


vis 


0.02 


+0.05886 


1 


52008.85 


2001.2700 


2 


ccd 


0.01 


-0.00388 


1 


52045.74 


2001.3710 


2 


ccd 


0.01 


+0.01421 


1 


52066.6993 .... 


2001.4283 


1 


ccd 


0.0008 


-0.00017 


1 


52432.88040 . . . 


2002.4309 


2 


ccd 


0.00050 


-0.00042 


4* 


52490.72613 . . . 


2002.5893 


1 


ccd 


0.00020 


-0.00027 


4 


52785.7012 .... 


2003.3969 


1 


ccd 


0.0010 


-0.00046 


5 


53154.4210 .... 


2004.4064 


1 


ccd 


0.0017 


+0.00027 


6 


53154.4212 .... 


2004.4064 


1 


ccd 


0.0015 


+0.00047 


6 


53209.7288 .... 


2004.5578 


1 


ccd 


0.0004 


+0.00021 


7 


53430.95992 . . . 


2005.1635 


1 


ccd 


0.00016 


-0.00011 


8 


53870.8862 .... 


2006.3679 


2 


ccd 


0.0005 


+0.00101 


9* 


53907.7573 .... 


2006.4689 


2 


ccd 


0.0002 


+0.00020 


9 


53928.7308 .... 


2006.5263 


1 


ccd 


0.0003 


+0.00002 


9 


54297.44944 . . . 


2007.5358 


1 


ccd 


0.00040 


-0.00041 


10 


54297.4498 .... 


2007.5358 


1 


ccd 


0.0008 


-0.00005 


11 


54297.4509 .... 


2007.5358 


1 


ccd 


0.0004 


+0.00105 


12 


54331.7839 .... 


2007.6298 


2 


ccd 


0.0003 


-0.00013 


13 


54368.6555 .... 


2007.7307 


2 


ccd 


0.0003 


-0.00043 


13 


54589.8873 .... 


2008.3364 


2 


ccd 


0.0002 


-0.00008 


13 


54647.7328 .... 


2008.4948 


1 


ccd 


0.0006 


-0.00016 


13 


54647.7332 .... 


2008.4948 


1 


ccd 


0.0002 


+0.00024 


13 


54868.9645 .... 


2009.1005 


1 


ccd 


0.0007 


+0.00010 


13 


Note. — Type: 1 — primary, 2 — secondary. Instrument: pg — photographic, vis — visual, 
ccd — CCD. References: 1 - ITorres ct al. (2001j[): 2 - Variable Star and Exoplanet Section of 



the Czech Astrono mica l Society tht tp : //var2 . astro . cz/EK/brno/ eclipsing_binaries . php) ; 3 - 
IZessewitschl (119541) ; 4 -ILacvl (2002): 5 - ILacvl }2003 | V; 6 -IHubs cher (200 5j); 7 - |Lacy| 1200411 : 
8 - ILacvl d2006fl : 9 - ILacvl rt2007T ): 10 - IBrat et al.l ll2007f) : 11 - IHubscher et alj d2008l) : 12 - 

Dicthclm (2008); 13 - This paper; * - Rcmcasurcmcnt. 
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Fig. 4. — Residuals from our fit to the primary and secondary 
eclipse timings of LV Her given in Table [4] 



Mount Hopkins, Arizona. We collected 42 spec tra with 
a photon-counting intensified Reticon detector (|Lathaml 
1985, 1992) in a single order 45 A wide centered near 
5187 A. The main features in this spectral window are 
the lines of the Mg I b triplet. The resolving power 
provided by this setup is A/AA » 35,000. Two addi- 
tional observations, for a total of 44, were obtained with 
a nearly identical system on the 1.5m Wyeth reflector at 
the Oak Ridge Observatory (Harvard, Massachusetts). 
The signal-to-noise ratios of these spectra range from 13 
to 28 per resolution element of 8.5 kms -1 . 

Radial velocities for the two components were derived 
using TO DCOR, a two-dimensi onal cross-correlation 
technique (jZucker fc Mazchj I1994J ) . This method uses 
two templates, one for each component of the binary, 
which we selected from a large library of synthetic spec- 
tra based on mode l atmospheres by R. L. Kurucz (see 
lLatham et al.l [20021 ). These templates have been calcu- 
lated for a wide range of effective temperatures (T e gp), 
surface gravities (\ogg), rotational velocities (vsini when 
seen in projection) , and metallicities ([m/H]). Following 
iTorres et al.l (|2002f ) the optimum templates for each star 
were determined by means of extensive grids of cross- 
correlations with TODCOR, seeking to maximize the 
average correlation weighted by the strength of each ex- 
posure. Because of the strong correlation between tem- 
perature, surface gravity, and metallicity in our narrow 
spectral window, we initially assumed solar metallicity, 
and surface gravities of logg = 4.0 for both components, 
close to the values determined from a preliminary anal- 
ysis. The projected rotational velocities we obtained are 
v sin i — 13 ± 1 kms -1 for both stars, and the best-fit 
temperatures were near solar, with a difference of only 
20-30 K between the primary and secondary. For the 
radial velocity measurements we adopted T e g = 5750 K. 
In §2] we describe experiments in which we extended the 
grids of correlations to other compositions in order to 
fine-tune the temperatures and attempt to constrain the 
metallicity. These minor changes in the stellar parame- 
ters have no effect on the radial velocities. Typical er- 
rors for our measurements are 0.8 kms -1 for the pri- 
mary and 1.1 kms -1 for the secondary. The stability of 
the zero-point of the CfA velocity system was monitored 
by means of exposures of the dusk and dawn sky, and 



TABLE 5 

Spectroscopic orbital solution for LV Her. 

Parameter Value 

Adjusted quantities 

P (days) a 18.4359535 

Ti (HJD-2,400,000) a 53,652.19147 

K x (kms" 1 ) 67.24 ± 0.19 

K 2 (kms" 1 ) 68.59 ± 0.27 

7 (kms" 1 ) -10.278 ± 0.094 

e 0.61273 ± 0.00073 

u) (dog) 352.20 ± 0.24 

Derived quantities 

Mi sin 3 i (M Q ) 1.193 ± 0.010 

M 2 sin 3 i (M Q ) 1.1697 ± 0.0080 

q = M 2 /Mi 0.9803 ± 0.0047 

ai sini (10 6 km) 13.471 ± 0.037 

a 2 sini (10 6 km) 13.743 ± 0.054 

asini (Rq) 39.119 ± 0.095 

Other quantities pertaining to the fit 

N oba 44 

Time span (days) 4018.0 

cri (kms" 1 ) 0.76 

<t 2 (kms- 1 ) 1.10 

a Ephcmcris adopted from § 12.21 

small systematic run-to- run correction s were applied in 
the manner described bv lLathaml (|1992t ). We determined 
also the light ra t io dir ectly from our spectra following 
iZucker fc M azehl (fl99l . and obtained i 2 /k = 0.95±0.03 
at the mean wavelength of our observations, which is es- 
sentially also the ratio in the visual band given that the 
stars are nearly identical in temperature. 
Although TODCOR significantly reduces systematic 
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Fig. 5. — Radial- velocity measurements for LV Her (filled circles 
for the primary, open for the secondary) along with the curves 
computed from our orbital solution. Phase 0.0 corresponds to the 
time of primary eclipse. The dotted line represents the velocity of 
the center of mass, and the error bars are smaller than the symbol 
size. O — C residuals for the primary and secondary are shown in 
the bottom panels. 
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TABLE 6 

Radial velocity measurements of LV Her. 



HJD 




RV 2 




°"2 


(O-C)i 


(o-c) 2 




(2,400,000- 


-) (kins- 1 ) 


(kms" 1 ) 


(kms -1 ) 


(kms -1 ) 


(kms 1 ) 


(km s _1 ) 


Phase 


48407.8794 . 


-36.77 


+16.94 


0.75 


1.09 


-0.38 


+0.59 


0.5389 


48428.8481 . 


-30.28 


+10.75 


0.77 


1.12 


+0.32 


+0.30 


0.6763 


48435.8126 . 


+2.41 


-24.64 


0.81 


1.17 


-0.56 


-0.85 


0.0541 


48459.7950 . 


-33.43 


+15.58 


0.84 


1.22 


+0.97 


+1.25 


0.3549 


48461.7375 . 


-35.70 


+18.08 


0.85 


1.23 


+0.88 


+1.53 


0.4603 


48462.6830 . 


-36.20 


+17.06 


0.84 


1.22 


+0.44 


+0.44 


0.5116 


51709.7370 . 


-33.05 


+12.51 


0.82 


1.19 


+0.18 


-0.62 


0.6377 


51710.7867 . 


-28.71 


+8.12 


0.69 


1.00 


+0.21 


-0.61 


0.6947 


51711.7546 . 


-22.10 


+0.39 


0.72 


1.04 


-0.34 


-1.05 


0.7472 


51712.7284 . 


-6.38 


-13.66 


0.81 


1.18 


+1.76 


-1.20 


0.8000 


51715.7770 . 


+62.12 


-84.03 


0.77 


1.11 


-0.20 


+0.31 


0.9653 


51740.7610 . 


-34.04 


+12.50 


0.64 


0.93 


-0.98 


-0.46 


0.3205 


51741.7446 . 


-35.18 


+14.19 


0.72 


1.04 


-0.18 


-0.75 


0.3739 


51742.7170 . 


-36.06 


+15.62 


0.74 


1.07 


+0.12 


-0.52 


0.4266 


51744.7057 . 


-37.70 


+15.90 


0.83 


1.21 


-1.26 


-0.51 


0.5345 


51799.6219 . 


-36.83 


+17.30 


0.73 


1.05 


-0.20 


+0.69 


0.5132 


51800.6010 . 


-37.55 


+14.68 


0.69 


1.00 


-1.65 


-1.17 


0.5663 


51802.6229 . 


-30.48 


+10.83 


0.66 


0.95 


+0.14 


+0.36 


0.6760 


51803.6558 . 


-24.42 


+3.48 


0.64 


0.92 


-0.15 


-0.51 


0.7320 


51831.5984 . 


-28.65 


+8.50 


0.75 


1.08 


+0.24 


-0.21 


0.2477 


51833.6157 . 


-35.73 


+14.79 


0.75 


1.09 


-1.25 


+0.38 


0.3571 


51834.5731 . 


-35.52 


+17.04 


0.72 


1.04 


+0.34 


+1.22 


0.4091 


51858.5764 . 


-26.02 


+5.88 


0.73 


1.06 


+1.09 


-1.01 


0.7110 


51861.5870 . 


+42.50 


-66.00 


0.77 


1.12 


-1.39 


-0.46 


0.8743 


51973.9772 . 


+56.88 


-78.61 


0.79 


1.14 


+0.26 


-0.09 


0.9706 


51977.9683 . 


-23.19 


+4.08 


0.68 


0.98 


+0.30 


+0.88 


0.1871 


52006.9280 . 


-20.18 


-0.38 


0.83 


1.21 


-0.49 


+0.30 


0.7579 


52009.8622 . 


+96.55 


-119.15 


0.81 


1.18 


+0.79 


-0.70 


0.9171 


52010.8700 . 


+54.74 


-77.11 


0.80 


1.16 


-0.68 


+0.19 


0.9717 


52011.9261 . 


+13.95 


-33.43 


0.76 


1.10 


+0.58 


+0.97 


0.0290 


52032.9120 . 


-21.81 


+0.64 


0.75 


1.09 


-0.63 


-0.20 


0.1673 


52034.8450 . 


-28.93 


+9.98 


0.71 


1.04 


+1.59 


-0.39 


0.2722 


52035.8449 . 


-33.46 


+14.00 


0.70 


1.02 


-0.15 


+0.78 


0.3264 


52039.9325 . 


-36.23 


+16.20 


0.73 


1.05 


+0.02 


-0.01 


0.5481 


52102.6290 . 


+80.99 


-102.64 


0.87 


1.26 


-0.07 


+0.81 


0.9489 


52138.6401 . 


+82.46 


-106.28 


1.31 


1.91 


+1.04 


-2.45 


0.9022 


52156.6725 . 


+51.60 


-73.90 


0.71 


1.04 


+0.17 


-0.67 


0.8803 


52157.6303 . 


+95.25 


-118.44 


0.70 


1.01 


-0.44 


-0.06 


0.9323 


52211.5883 . 


+27.39 


-42.70 


0.89 


1.30 


+0.01 


+5.99 


0.8590 


52360.9241 . 


+69.40 


-91.63 


0.78 


1.13 


+0.24 


-0.31 


0.9593 


52420.8295 . 


-25.10 


+4.20 


0.77 


1.11 


+0.58 


-1.23 


0.2087 


52421.8579 . 


-30.47 


+9.39 


0.78 


1.13 


-0.43 


-0.49 


0.2645 


52423.9196 . 


-35.64 


+15.54 


0.73 


1.06 


-0.57 


+0.53 


0.3763 


52425.8862 . 


-35.66 


+17.17 


0.77 


1.11 


+1.03 


+0.51 


0.4830 



Note. 
text. 



Radial velocities arc in the heliocentric frame, and include all corrections described in the 



errors in the radial velocities caused by line blending, 
residual effects can remain as a result of shifts of the 
spectral lines in and out of our narrow spectral window 
as a function of orbital phase. We investigated these ef- 
fects by mea ns of numerica l simul ations si milar to those 
described by [Latham et al. (199(f) (see also iTorres et al.l 
119971 l2000al ). We generated synthetic composite spectra 
matching our observations by combining copies of the pri- 
mary and secondary templates used above, shifted to the 
appropriate velocities for each of the exposures as pre- 
dicted by a preliminary orbital solution, and scaled to the 
observed light ratio. These synthetic observations were 
then processed with TODCOR in exactly the same way 
as the real spectra, and the resulting velocities were com- 
pared with the input shifts. The differences for LV Her 
were typically well under 0.5 kms -1 for both stars, but 
were nevertheless applied as corrections to the raw veloc- 
ities. They affect the masses at the level of 1.3% for the 
primary and 0.8% for the secondary, which are similar 
to the formal errors. The corrected measurements are 



listed in Table [5] along with their uncertainties. Similar 
adjustments based on the same simulations were made 
to the light ratio, and are already included in the value 
reported above. 

Our orbital solution is shown in Figure [5l and was de- 
rived holding the ephemeris fixed according to § 12.21 The 
residuals are shown as well, and are listed in Table [Sj 
The orbital elements and derived quantities (minimum 
masses, semimajor axes, etc.) are given in Table [SJ 

3. MODELING OF THE PHOTOMETRIC OBSERVATIONS 

In order to remove the small nightly offsets in the dif- 
ferential photometry, we performed preliminary fits of 
the URSA and NF O light curves by using the Nelson- 
Davis-Etzel model (popper fc Etzeii [MlTlEtzel HMD) 
as implemented in the EBOP code. This model is 
well suited for well-detached systems such as LV Her. 
The nightly residuals, which were typically less than 
0.02 mag and were uncorrelated with phase, were re- 
moved from each light curve. This improved the resid- 
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TABLE 7 

Photometric orbital solutions for. LV Her. 



Parameter URSA NFO Adopted 



J 2 0.967 ± 0.004 0.983 ± 0.003 0.9772 ± 0.0080 

n+r 2 0.06880 ± 0.00009 0.06806 ± 0.00006 0.06829 ± 0.00037 

k = r 2 /n 0.986 ± 0.015 0.961 ± 0.009 0.968 ± 0.012 

n 0.03465 ± 0.00026 0.03471 ± 0.00015 0.03470 ± 0.00028 

r 2 0.03415 ± 0.00026 0.03335 ± 0.00016 0.03357 ± 0.00028 

i (deg) 89.489 ± 0.010 89.509 ± 0.009 89.500 ± 0.010 

ecosu 0.60662 ± 0.00011 0.60739 ± 0.00006 0.60721 ± 0.00038 

esinw -0.0938 ± 0.0019 -0.0796 ± 0.0012 -0.0836 ± 0.0071 

e 0.61383 ± 0.00018 0.61259 ± 0.00010 0.61288 ± 0.00062 

a; (deg) 351.21 ± 0.18 352.53 ± 0.12 352.12 ± 0.66 

h 0.516 ± 0.007 0.524 ± 0.004 0.5220 ± 0.0040 

£ 2 0.484 ± 0.007 0.476 ± 0.004 0.4780 ± 0.0040 

(£ 2 /h)v 0.938 ± 0.024 0.908 ± 0.014 0.916 ± 0.015 

u 1 = u 2 0.53 ± 0.04 0.55 ± 0.04 0.54 ± 0.03 

N obs 6690 2946 

Time span (days) .... 2970.9 1514.9 

a v (mmag) 12.49 6.42 



uals of the URSA data by about 27%, and those from 
NFO by a more significant 42%. The corrected light 
cu rves were then fit t ed by using the JKTEBOP code 
of iSouthworth et all (|2007h . based on the same model, 
which allows for more realistic estimates of the uncer- 
tainties. The adjustable parameters are the central sur- 
face brightness J 2 of the secondary relative to the pri- 
mary, the inclination angle i, the sum of the relative 
radii r\ + r 2 , the ratio of the radii k = r 2 /r± , the eccen- 
tricity factors ecosw and esinw, a phase shift, and the 
magnitude at quadrature. A linear limb-darkening law 
was adopted, consistent with our experience that with 
the amount and precision of our data, nonlinear laws do 
not improve the accur acy of the fits significantly (e.g., 
lLacv et al.l 120051 I2008D . The coefficient u was left free, 
and constrained to be the same for the two stars. The 
gravity brightening exponent was s et to 0. 3 5 for both 
stars, based on the calculations by IClaretl (|1998f ) and 
the mean temperature and surface gravity of the com- 
ponents. The mass ratio q was adopted from the spec- 
troscopy. 

Separate solutions for the URSA and NFO photom- 
etry are given in Table [7J The uncertainties for the 
fitted parameters were evaluated by Monte Carlo sim- 
ulation with JKTEBOP. Tests indicated third light was 
not significant, and it was therefore set to zero. The 
last column of the table gives the adopted fit (weighted 
average), with uncertainties that account for the differ- 
ence between URSA and NFO. The uncertainties for the 
individual radii were derived followin g the prescription 
recommended bv lTorres et al.l (|2000bD . based on the er- 
rors for r\ + r 2 and k. The synthetic curve corresponding 
to the adopted solution is displayed in Figures [T]-[3] along 
with the observations. The eccentricity and u> from the 
light curves are in excellent agreement with the spectro- 
scopic values in Table [51 and the photometric light ratio 
{(■2/(-i)v = 0.916 ± 0.015 is consistent with the estimate 
of £ 2 /£ 1 = 0.95 ± 0.03 from the CfA spectra. The well- 
detached stars are essentially spherical, and both eclipses 
are partial, with 88.4% of the primary light blocked at 
the primary minimum, and 95.6% of the secondary light 
covered at the secondary minimum. 

We note that the linear limb-darkening coefficient we 
determined, u — 0.54 ± 0.03, is consistent with theo- 



retical value of 0.5 68 for the visual band according to 
Ivan Hammd (|1993[ ). but is signi ficantly small er than the 
value 0.658 from the tables bv iClareti ((2000), for a star 
of this temperature and surface gravity. Similar discrep- 
ancie s have been found in other systems (see, e.g.- IClaretl 
2008). As a test we repeated the light curve fits holding 
the limb-darkening coefficient fixed at this last theoreti- 
cal value. The change in r\ + r 2 was only +0.15%, the 
individual radii increased by less than 0.1%, and the in- 
clination angle decreased by 0?065. These changes do 
not impact the absolute masses or radii at a significant 
level, compared to other sources of uncertainty. 

4. ABSOLUTE DIMENSIONS 

The photometric orbit and the spectroscopic orbit may 
be combined to yield the absolute masses and radii of the 
binary stars' components. These are given in Table [51 
along with other properties described below. The masses 
and radii are determined to 0.9% or better for both stars. 
The tests and checks for systematics described above sug- 
gest that they are also accurate at this level. 

The ratio of the effective temperatures is very well con- 
strained from the light curves through the central surface 
brightness parameter J 2l a ccording to w hich the stars 
differ by only about 30 K (|Popperl [l980l Table 1). A 
similar difference was obtained from the CfA spectra. 
The absolute temperatures of the stars may in principle 
be estimated also from the spectroscopy, though in prac- 
tice this is difficult with the material at hand because of 
the correlation with metallicity mentioned earlier. We 
first determined the magnitude of this dependence by 
repeating the grids of correlations described in § 12.31 for 
metallicities between [m/H] = —1.0 and [m/H] = +0.5 
in steps of 0.5 dex, and surface gravities of log g = 4.0 
and 4.5. At each composition we determined the primary 
and secondary temperatures interpolated to the surface 
gravities indicated in Table [8j as well as the luminosity- 
weighted mean temperature for the system (where the 
weights depend only on the radius ratio and temperature 
ratio that come from the light curve fits). An additional 
constraint on the mean system temperature is available 
from the absolute photometry of LV Her presented ear- 
lier, through color indices and color-temperature calibra- 
tions (including metallicity terms). Equating the two 
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TABLE 8 
Physical properties of LV Her. 



Parameter Primary Secondary 



Mass (M ) 


1.193 


± 0.010 


1.1698 


± 


0.0081 


Radius (Rq) 


1.358 


± 0.012 


1.313 


± 


0.011 


log g (cgs) 


4.2493 


± 0.0082 


4.2695 


± 


0.0081 


Temperature (K) 


6060 


± 150 


6030 


± 


150 


logL (L Q ) 


0.349 


± 0.044 


0.311 


± 


0.044 


M bol (mag) 

\ ' 


3.86 


± 0.11 


3.95 


± 


0.11 


-0.04 


± 0.10 


-0.04 


± 


0.10 


My (mag) b 


3.90 


± 0.16 


4.00 


± 


0.16 


a (R @ ) 




39.120 ± 0.095 






Distance (pc) 




352 ± 24 






[m/H] 




+0.08 ± 0.21 






vs'mi (kms - 1 ) c 


13 


± 1 


13 


± 


1 


fpsyncsini (kms _1 ) d ... 


16.0 


± 0.1 


15.5 


± 


0.1 


"peri sin i (kms - 1 )° 


19.6 


± 0.2 


19.0 


± 


0.2 



a BCv is taken from Flower (1996), and the uncertainty includes the 
contribution from the temperature as well as an additional 0.10 mag 
added in quadrature. 

The bolometric magnitude adopted for the Sun is M° ol = 4.732, 
for consistency with the zero-point of the bolometric corrections from 
IFlowerl U996T) . 

c Value measured spcctroscopically. 

^ Projected pseudo-synchronous rotational velocity. 

° Projected synchronous rotational velocity at periastron. 

mean temperatures then allows to solve for the metallic- 
ity. With the photometry listed in Table [3] we formed 7 
different color indices, we de-reddened them as indicated 
in 5 12. f I an d we made u se of t h ree different calibra t ions b y 
Rami rez k, Melendezl (12005ft. iCasagrande et all (|2006f ). 
and lGonzalez Hernandez fe Bonifacio! (|2009f h While the 
7 indices are not independent, they are useful in pro- 
viding an idea of the internal consistency of the results. 
The agreement between the different indices and differ- 
ent calibrations is excellent, generally within 100 K. The 
metallicity dependence of the photometric temperatures 
is small, but was nevertheless accounted for by itera- 
tions with the spectroscopic results. We obtained a mean 
system temperature of T c g = 6050 ± 140 K and [m/H] 
= +0.08 ± 0.21, where the errors are dominated by the 
uncertain reddening. The primary and secondary tem- 
peratures are then 6060 ± 150 K and 6030 ± 150 K, re- 
spectively, though of course AT c ff is known much more 
precisely. These temperatures correspond to a spectral 
type of approximately F9 (e.g.. lPopper| [i~980). 

L V Her does not have a n entry in the Hipparcos cata- 
log (jPerrvman et al.l 119971 ). The distance to the system 
is estimated here as 352 ± 24 pc, based on the bolo- 
metric luminosities, the apparent system magnitude of 
V = 11.02 ± 0.02 (Table EJ), the adopted extin ction of 
3.1 x E (B — V), and bolometric corrections from IFlowerl 
(1996). Separate distance estimates for each component 
agree nearly perfectly, showing the internal consistency 
of the results. 

Also listed in Table [5] are the rotational velocities (pro- 
jected along the line of s ight) expe cted if the stars were 
pseudo-synchronized (see lHutJ[l981f ). as well as those that 
correspond to synchronization with the angular velocity 
at periastron, for comparison with the measured vsmi 
values. We discuss these below. 

5. COMPARISON WITH THEORY 

The high accuracy and precision of our measurements 
for LV Her offers the opportunity to test various aspects 



of theoretical modeling. In this particular case, knowl- 
edge of the chemical composition, even though the pre- 
cision is not as high, is especially interesting since it is 
rarely available for eclipsing binaries and it eliminates 
one of the free parameters in the comparison with the- 
ory. 

5.1. Stellar evolution theory 

In Figure we compare the two best determined pa- 
rameters of LV Her (M and R) against mode l isochrones 
from the Yonsei-Yale series bv lYi et all ()2001[ ) , computed 
for the metallicity of [Fe/H] = +0.08 determined in §SJ 
These models include convective core overshooting (see 
iDemarque et al.ll2004f) . which has an effect for stars of 
this mass and higher, and treats convection in the stan- 
dard mixing length approximation, with a mixing length 
parameter «ml = 1-7432, calibrated against the Sun. 
An isochrone corresponding to 3.85 Gyr (solid line) pro- 
vides an excellent fit, and the high precision of the mea- 
surements suggests an age uncertainty of no more than 
0.4 Gyr for this fixed metallicity. However, the rela- 
tively large error in the measured [Fe/H] adds signif- 
icantly to the age uncertainty. This is shown by the 
shaded area around the 3.85 Gyr isochrone corresponding 
to o~ [Fe/H] — 0-21 dex, which contributes at least another 
0.5 Gyr. In the lower panel of the figure the temperature 
is shown as a function of mass, with the same models as 
in the top panel. Once again the isochrone of 3.85 Gyr 
matches the observations very well, but in this case the 
lower precision of T c g provides a much weaker constraint 
on the age, and gives rise to a larger uncertainty than 
the error that comes from [Fe/H]. 

The T ff-rog g diagram in Figure [7] implicitly compares 
all four measured quantities (M, R, T c s, and [Fc/H]) 
against evolutionary tracks from the Yonsei-Yale series. 
The models are computed for the observed metallicity 
of [Fe/H] = +0.08, and for the precise masses we mea- 
sure. The shaded area indicates the uncertainty in the 
location of the tracks that comes from the mass errors. 
The same 3.85 Gyr isochrone as before is shown with a 
dashed line. The stars are seen to be in the middle of 
the main-sequence band. 

A similarly good fit to the properties of LV He r is prq - 
vided by the stellar evolution models of Claret (2004), 
shown in Figure also for a metallicity constrained to 
the measured value. These models are similar to the 
ones used before in their treatment of convection (oml = 
1.68), and also include overshooting (a ov = 0.20), but 
differ in the details. In this case, we find that the best-fit 
age is 4.2 Gyr, about 9% older than with the Yonsei-Yale 
models. The difference is of the same order as the age 
uncertainty discussed above. 

5.2. Tidal theory 

The high orbital eccentricity and our accurate mea- 
surement of the projected rotational velocities of both 
components constitute useful probes of tidal forces in 
the system and allow for an interesting comparison with 
the predictions from theory. For our initial compa risons 
we ha ve us ed the radiat ive damping formalism of iZahnl 
(| 19771) and IZahnl ([1989), often invoked for this type of 
analysis, to calculate the critical times for LV Her corre- 
sponding to synchronization of the axial rotations, align- 
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merit of the rotational and orbital axes, and circular- 
ization of the orbit. These are then compared with the 
evolut ionary age, w hich is 4.2 Gyr according to the mod- 
els by iClaretl 120041 used in this section. The procedure 




£ 5500 - 



Fig. 6. — Measurements for LV Her compared aga inst the Yonsei- 
Yale models flVi et al.H200HIDemarque et al.ll2004l) . (a) Radius vs. 
mass, shown with isochrones from 1 Gyr to 6 Gyr in steps of 1 Gyr 
(dotted lines) for a fixed metallicity of [Fe/H] = +0.08. The best- 
fit isochrone for 3.85 Gyr is indicated with a solid line. The shaded 
area indicates the uncertainty that comes from the 0.21 dex error 
in the measured metallicity, at this particular age. (b) Effective 
temperature vs. mass. 
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Fig. 7. — Surface gravity vs. effective temperature diagram for 
LV Her. The solid l ines correspond to evolutionary tracks from the 
Yonsei-Yale series HYi et al.l f2001: Dcmaraue ct al. 200D for the 
exact masses we measure, and for the observed chemical composi- 
tion. The effect of the mass uncertainties in the location of these 
tracks is indicated by the shaded areas. The dashed line represents 
the best-fit isochrone for 3.85 Gyr. 
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Fig. 8. — Surface gravit y vs. ef fective temperature of LV Her 
compared with models by Claret (2003). The solid lines corre- 
spond to evolutionary tracks for the measured masses, and for the 
observed chemical composition. The age according to these models 
is 4.2 Gyr. 



follows closely that described bv I Claret &: Cunhal <|1997h . 
and consists of integrating the differential equations de- 
scribing the changes in these quantities along the evolu- 
tionary track for each star. 

We find that rotational synchronization due to tidal 
forces is not expected to occur for the primary un- 
til the system reaches an age of r synCj i = 5.7 Gyr 
(logTsynd = 9.754), and the secondary synchronization 
is predicted to happen shortly after (T synCi 2 = 5.9 Gyr, 
logT synCi 2 = 9.774). This appears to be in agreement 
with the fact that the measured v sin i values are slower 
than expected for pseudo-synchronization in the eccen- 
tric orbit (Table [8]), assuming the spin axes are paral- 
lel to the orbital axis (but see below). The time of 
orbital circularization indicated by theory is 7.3 Gyr 
(logT C i rc = 9.862), again consistent with the fact that 
the orbit is observed to be highly eccentric. These times 
are displayed graphically in Figure [9l along with the evo- 
lution of the radius of each star according to the evolu- 
tionary tracks. 




9.7 

Log age (yr) 

Fig. 9. — Radius as a function of age for LV Her (solid line for 
the prima ry, dashed for the secondary) from the stellar evolution 
models by Claret ( 2004]) , for the measured masses and metallicity. 
Measured radii and their uncertainties are represented by the hor- 
izo ntal s h aded bands . The critical times according to the theory 
by Zahn (1977, 1989) are indicated with arrows for synchroniza- 
tion and spin-orbit alignment of each component, as well as for 
circularization of the orbit. The evolutionary age is also shown. 
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A somewhat unexpected result is that the times of 
alignment predicted by theory are also significantly 
longer than the current age of the system: r a i; gni i = 
6.4 Gyr and T a ii g n,2 = 6.8 Gyr (or logr a ii g n,i = 9.804 and 
logTaiign,! = 9.832, in reference to Figure [5]). At face 
value this indicates both spin axes may be inclined rela- 
tive to the axis of the orbit, which is the opposite of what 
is almost universally assumed for binary systems. The 
long orbital period of LV Her could provide a plausible 
explanation, in principle, but in any case the misalign- 
ment makes it somewhat problematic to interpret the 
apparent lack of (pseudo-)synchronization inferred from 
the measured usini, since the projection factor can no 
longer be assumed to be known. 

The calculations above involve a number of approxi- 
mations implicit in th e equa tions we have us ed, as de- 
scribed, e.g., bv iZahnl l| 19771 ) and [Hut] (|198lT ). In par- 
ticular, they are linearized around the equilibrium state, 
and are strictly valid only for relatively small eccentric- 
ities and near-synchronous rotation, and for small rela- 
tive inclinations between the spin axes and the axis of 
the orbit. We have also ignored changes in the semi- 
major axis that occur concurrently with the evolution of 
other orbital elements. The condition on the eccentricity 
is most certainly not met for LV Her (e ~ 0.613), and 
it is unclear a priori to what extent this may affect the 
conclusions above. This, and the potentially interesting 
situation regarding the spin axes, have prompted us to 
explore the situation in greater detail. We have directly 
integrated the more general differential equations (valid 
also for high eccentricities) that describe the evolution of 
the semimajor axis (da/dt), eccentricity (de/dt), angu- 
lar rotational rates (dCli/dt, d^/dt), and inclination of 
the spin axes relative to th e orb i tal ax is (dAi/dt in our 
nomenclature) as given by iHutl (|1981l ). using a fourth- 
order Runge-Kutta method. Because they are coupled, 
these equations must be integrated simultaneously, and 
as before the stellar properties have been interpolated 
from the evolutionary tracks at each time step. The tur- 
bulent dissipation timescale for stellar phases with con- 
vective envelopes is taken to be (MR^/L) 1 ' 3 , where M, 
R, and L are the mass, radius, and luminosity of the 
star. For phases with radiative envelopes the timescales 
adopted follow closely those in eq.(17) and eq.(18) by 
IClaret fc Cunhal $39% . 

The initial conditions are of course not known, so we 
explored a range of values that give predictions matching 
the observed quantities at the current age of the system. 
This is shown in Figure HOa .b for the eccentricity and 
orbital period. The initial values are eo = 0.66 and Pq = 
21.8 days at the starting age of log To = 7.0. According 
to these calculations we expect the orbit to circularize at 
an age logr C i rc ~ 9.85 (r C j rc ~ 7.1 Gyr), which is in fairly 
good agreement with our previous estimate despite the 
concerns expressed above. 

Also shown in the figure is the evolution of the angular 
rotation rate Q = 27r/P ro t for each star, which we nor- 
malize for convenience to the orbital rate £l rb — 27r/P or b 
(FigurefTOb.d; solid lines). For LV Her we have no obser- 
vational constraints on fii and f^2, so the initial values 
are arbitrary and chosen to be the same. In principle they 
could be constrained by directly measuring the rotation 
periods of the stars due, e.g., to rotational modulation by 



spots, although in practice this is difficult for a spatially 
unresolved binary. More importantly, these stars may be 
too hot to be affected significantly by spots. The ratio 
between the pseudo-synchronous velocity and the mean 
orbital motion is shown in Figure 110b . and this curve 
is repeated in Figure flOb.d (dot-dashed line). As seen 
from the convergence of the dot-dashed and solid curves, 
the stars are predicted to reach pseudo-synchronization 
roughly at logr psync ~ 8.75 (r psync ~ 0.6 Gyr), which is 
much younger than the evolutionary age. This is com- 
pletely at odds with our previous estimate, showing the 
limitations of that approach for a system like LV Her. 

We also have no direct constraint on the relative angle 
Ai between the spin axes and the orbital axis. Further- 
more, the relation between Ai and the orbital and rota- 
tional inclinations i or b and i ro t, both measured relative 
to the plane of the sky, is given by 

cos Ai = cos iorb cos £ ro t + sin i orD sin i rot cos A , (1) 

which involves an unknown angle A between the sky- 
projected angular momentum vectors of the orbit and 
the stellar spin. 

The spectroscopically measured projected rotational 
velocities of the stars, which we refer to more properly 
now as v sin i rot , do provide some indirect constraint on 
a combination of the theoretically predictable quantities, 
but this still involves the unknown angle A. Because i or b 
is very nearly 90° for LV Her (see Table [S]), we may make 
the approximation cosAi w sini rot cos A. We then have 



«l 2 sini ro t 



2tt cos Ai 

-Porb ^orb COS A 



R 



(2) 



In this equation all quantities on the right-hand side are 
either known from stellar evolution calculations (R), or 
can be computed from the solution of the differential 
equations, except for cos A, which depends on the ob- 
server's viewpoint. For the sake of illustration, we ignore 
this term in the following (or consider A to be small), so 
that cos Ai w sin i rot . In Figure fTOF we show the evolu- 
tion of Ai for four different initial values (20°, 40°, 60°, 
80°). Curves for the primary and secondary are nearly 
indistinguishable, so we show only those for the primary. 
Once again the conclusion is very different from the one 
we reached with the equations by Zahn: alignment of 
the spin axes with the orbital axis is achieved much ear- 
lier than the age of the system, at logT a n gri ~ 8.9 (or 
Taiign — 0.8 Gyr). The corresponding usini rot curves for 
the primary and secondary, calculated using eq.([2]), are 
shown in the bottom panels and are compared with the 
measured rotations (dashed line and shaded uncertainty 
region). The general trend is for the stars in LV Her to 
be spun up by tidal forces, but the calculations disagree 
with the observations at the current age of the system, 
predicting the measured rotations should be 16 kms -1 
for the primary and 15.5 kms -1 for the secondary (see 
also Table IS]). 4 These are nominally outside of the range 
allowed by the observational uncertainties. We consider 
it very unlikely that our measurements are off by as much 

4 Although the rotation rates Q± and O2 have been chosen arbi- 
trarily here, this does not affect the predicted rotational velocities 
at the current age because pseudo-synchronization occurs much 
earlier, independently of the starting values of Hi and Q2- 
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Fig. 10. — Observations of LV Her compared against tidal theory 
IIHutlll98lI) . (a) Eccentricity evolution. The dashed line represents 
the measured value, and the vertical dotted line in this and the 
other panels marks th e current evol utionary age of 4.2 Gyr, accord- 
ing to the models by Claret (2004). (b) Evolution of the orbital 
period, with the current value represented with a dashed line, (c) 
Rotation rate of the primary, normalized to the orbital rate. The 
dot-dashed line represents the evolution of the pseudo-synchronous 
rate shown in panel (e). (d) Same as (c), for the secondary, (e) 
Evolution of the pseudo-synchronous rotation rate, (f) Evolution 
of the relative angle Ai between the spin axis and the orbital axis, 
for different initial values, (g) Evolution of the projected rotational 
velocity of the primary, for the same four trial values of Ai. The 
dashed line represents the measured value, its uncertainty indicated 
by the shaded area, (h) Same as (g), for the secondary. 

as 3 kms -1 for the primary and 2.5 kms" 1 for the sec- 
ondary, since our formal uncertainty of ±1 kms -1 is 
already conservative. The alternative would be a defi- 
ciency in the calculations, related to remaining approx- 
imations in the differential equations for this very chal- 
lenging problem. A combination of both effects is also 
possible. An even more general treatment of the dynam- 
ical problem of tidal friction in binary syste ms with de- 
forma ble components has been developed bv I Alexander! 
(|1973f ). but the application is considerably more complex 
and is beyond the scope of the present work. 

We note, finally, that an empirical test of the prediction 
of alignment is available in eclipsing bina ries by measur- 
ing the Rossiter-McL aughlin effect (e.g.. lAlbrecht et al.l 
2007; Al bredul [20081) . This provides a direct measure of 
the angle A between the sky projections of the spin axes 
and the axis of the orbit. 

6. DISCUSSION AND CONCLUSIONS 

As a result of our intensive spectroscopic and pho- 
tometric monitoring, the absolute masses and radii of 



LV Her are now determined to 0.9% or better, and are 
among the best available for any eclipsing binary. In ad- 
dition we have established the effective temperatures, as 
well as the overall metallicity. The latter quantity has 
not been determined for many eclipsing binaries, and is 
important because it reduces the degrees of freedom in 
the comparison with theory. Unfortunately the precision 
of [m/H] (and T c g to some extent) is limited in this case 
by poor knowledge of interstellar extinction in the direc- 
tion of the star. This could perhaps be remedied with 
additional absolute photometric observations. Neverthe- 
less, the combination of M, R, T Q g, and [m/H] provides 
unusually strong contraints on stellar evolution theory. 
Comparison with current models yields an excellent fit 
to the observations, and indicates the stars are approx- 
imately half way through their main-sequence phase, at 
an age between 3.8 and 4.2 Gyr, depending on the model. 

Eclipsing binaries with periods as long as that of 
LV Her are relatively rare and difficult to find, except per- 
haps in the course of automated variability surveys such 
as those designed to search for transiting planets, which 
often operate almost continuously for weeks or months 
at a time. The very high eccentricity of the system 
(e ~ 0.613) makes it the most ex treme case among those 
with well-determined properties ((Torres e t al. 2009]). Al- 
though apsidal motion is expected, the apsidal period is 
predicted to be very long (U ~ 58,000 yr) on account 
of the wide separation between the stars; it may take 
decades for the effect to be detected and measured accu- 
rately. 

Tidal theory is consistent with the observations in pre- 
dicting that the orbit should not yet have been circular- 
ized by tidal forces. Direct integration of the cou pled 
differe ntial equations for this problem according to iHutl 
(1981) indicates that the spin axes have already been 
forced into alignment with the orbital axis early on in 
the evolution of this system. The spins of both stars are 
also expected to be pseudo-synchronized with the orbital 
motion, although this is not quite in agreement with the 
measured rotations, and may indicate either large errors 
of measurement (which we believe are unlikely) or short- 
comings in theory, or some combination of both. The 
case of LV Her highlights the importance of bearing in 
mind the assumptions under which the equations for tidal 
evolution have been derived. Because of their ease of use 
and s i mpler interpret ation, the l inearized equations by 
IZahnl (p7l (see also iHul tl98ll IClaret fc CuriLl [199% 
have occasionally been applied to binary systems that 
may violate those assumptions to some degree (namely, 
the requirement of small eccentricities and spin-orbit 
inclinations, and near-synchronous rotations). Exactly 
when and how the theory breaks down is difficult to pre- 
dict. With its high eccentricity, LV Her is one such case, 
as we have demonstrated here by comparing the critical 
times from the simpler formulation with the res ults from 
the more general equations given bvlHutl (Il98l . 
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